FOCUSING OF SPHERICAL NONLINEAR PULSES IN M 1+3 , III. 
SUB AND SUPERCRITICAL CASES 

REMI CARLES AND JEFFREY RAUCH 



Abstract. We study the validity of geometric optics in L°° for nonlinear 
wave equations in three space dimensions whose solutions, pulse like, focus at 
a point. If the amplitude of the initial data is subcritical, then no nonlinear 
effect occurs at leading order. If the amplitude of the initial data is sufficiently 
big, strong nonlinear effects occur; we study the cases where the equation is 
either dissipative or accretive. When the equation is dissipative, pulses are 
absorbed before reaching the focal point. When the equation is accretive, the 
family of pulses becomes unbounded. 



1. Introduction 



This paper is the last of a series of three, after [CR02] and [CR04 . In these 
three papers, we consider the asymptotic behavior as s — > of solutions of the 
initial value problem 

(t,x) 6 [0,T] x M 3 , 



f Du e + a\d t u e \ p ~ 1 d t u e 



(1.1) 



0, 



«1t=o - z J+1U olr, 



r-r 



d t u e 



e J Ui r, 



£ 

r-r 



where □ := df — A x , a is a complex number, r = \x\, tq > 0, and, 1 < p < oo. The 
functions Uq and U\ are real-valued, infinitely differentiable, bounded, and there is 
a zq > so that for all r > 0, 



(1.2) 



SUpp£/j(r, .) C [-Zq, Zq]. 



The last assumption implies that at time t = the solutions are families of spherical 
pulses supported in a 0(e) neighborhood of r = r$. The initial data are spherically 
symmetric, thus in the limit e — > 0, a caustic is formed, reduced to the focal point 
(t, x) = (ro,0). Before going further into details, we rescale our parameters as in 
| CR.04 1 . Introduce e~ J u e —: u £ instead of u £ so that the solutions have derivatives 
of order Q(l ) away from the caustic. Define a := (p — 1) J. The initial value 



problem (1.1) is transformed to 



Du £ + ae a \dtu £ \ p " 1 d t u e = 0, (t,x) e [0,T] 



(1.3) 



lt=o 



= eU r 



r - r 



ro 



In | CR04 1 , formal arguments, inspired by the linear case a = 0, led to the following 
distinctions, in the spirit of those computed formally in | HK87 |, 
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a + 2 > p 


a + 2 = p 


a + 2 < p 


a > 


linear caustic, 
linear propagation 


nonlinear caustic, 
linear propagation 


supercritical caustic, 
linear propagation 


a = 


linear caustic, 
nonlinear propagation 


nonlinear caustic, 
nonlinear propagation 


supercritical caustic, 
nonlinear propagation 



In [CR02], we studied the case "linear caustic, nonlinear propagation" ; we proved 
that nonlinear geometric optics provides a good approximation of <9tu 6 away from 
the focal point (t, r) = (ro,0), and that the nonlinear term is negligible near the 
focus. In [CR04], we analyzed the case "nonlinear caustic, linear propagation". 



In some sense, it is the exact opposite of the previous case; the nonlinear term is 
negligible outside the focal point, but has a relevant influence near (r ,0), which 
is described by a nonlinear scattering operator. Moreover, this scattering operator 
broadens the pulses (at least if Uq and U\ are small), which leave the focus with 
algebraically decaying tails. 

In this paper, we discuss the remaining cases of the abov e ta ble. In the last 
three cases, we treat only the case of a real, that is when Eq. (1.3) is dissipative or 
accretive (in particular, we do not treat the case of conservative equations). 

The first case, "linear caustic, linear propagation" , suggests that the nonlinear 
term is everywhere negligible; we prove that this is so. In the last three cases, we 
assume that the coupling constant a is real. For the "supercritical caustic" cases, 
strong nonlinear effects are expected near the focal point. We prove that when 
the equation is dissipative (a > 0), then the dissipation is so strong near the focus 
that the pulses are absorbed. This result is the pulse analogue of JMR95 | and 
| JMR00 1 , which proved absorption in the case of wave trains (the initial profiles 



Uj are assumed to be periodic with respect to their last variable instead of com- 



pactly supported). More precisely, in ]JMR95fl and pMROCfl , it is proved that the 
exact solution u s of the dissipative (a > 0) wave equation ([Ll]) with J = 0, is 
approximated as follows, 



d t u e {t,x) ~ d t u(t, x) + U- \ t, x 



t+\x\ 



U 4 



t-\x\ 



where the profiles U± are periodic with respect to their last variable, with mean 
value zero. The absorption of oscillations is given by U± = past the caustic. Thus, 
only the average term remains, included in dtu. For an almost periodic function, 
the notion of average is given by 

i-T 



f 



lira — 

T^+oo 2T 



f(9)d9. 



—T 



When / is compactly supported, the above limit is zero, and pulses formally have 
mean value zero. Thus, the absorption of pulses is the formal analogue of the 
absorption of oscillations. 

Our present framework makes it possible to analyze very precisely the corre- 
sponding phenomenon for pulses; they are absorbed when approaching the caustic, 
that is even before reaching it. We prove that this phenomenon occurs for the last 
three cases of the table, "supercritical caustic, linear/nonlinear propagation" and 
"nonlinear caustic, nonlinear propagation". 

The present paper along with [CR02| and CR04[ prove that the distinctions 
derived formally in [CR04| and recalled in the above tables are correct. Let us 
give an interpretation of these results when the nonlinearity is fixed, and when one 
modulates the amplitude of the initial data in (1.1). Consider a fixed p > 2, and 
modify the value of J. For a unified, complete, presentation, we assume that the 
equation is dissipative, a > 0. 
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• If J > 2_£ j then the pulse is not affected by the nonlinearity at leading 
order. It remains too small to ignite the nonlinearity. 

• If J = > then the nonlinear term is negligible away from the focus, but 
the caustic crossing, described by a scattering operator, has enlarged the 
support of the pulse, and decreased its amplitude. The pulse is too small 
to see the nonlinearity outside the focal point, but the amplification near 
the caustic makes the nonlinear term relevant there. 

• If J < jEj, then the pulse is absorbed at the focus. It is sufficiently big to 
make the nonlinear effects so strong that the dissipation is complete before 
the focus. 

When p = 2, the nonlinear term is negligible if J > 0, and if J = 0, the pulses are 
absorbed. If 1 < p < 2, the same method would prov e that the nonlinear term is 
negligible if J > 0, and the case J = was treated in | CR02 |. 



B efore st ating precisely our results, we make a change of unknown, as in [CR02] 
and CR04 1 . Since the initial data are spherical, so is the solution. With the usual 
abuse of notation, 



u £ (t,x) = u £ (t,\x\), u e {t, \x\) e C, 
Introduce v £ := (uf_,u!L) where 

(1.4) u £ {t,r) :=ru £ {t,r), v% := {d t ± d r )u e , 



even m r \ 



r)- 



v £ e C°°(! 



Then (1.3) becomes 



(1.5) 



where 



r - r 



g(y) 



-a2- p \y\ p -\ 



T\t=0 



Px 



± ePi r, 



r - r 



P T (r,z) := rUi(r, z) ± rd z U (r, z) 



P\{r,z) := U (r, z) + rd r U (r, z) . 
We prove asymptotics for v £ ; asymptotics for dtu £ are deduced by (1.4), 

2r 

For the sub-critical case, introduce the solution of the linear equation 

(9 t ± 9r)(w±)froo = 0, 
((fi)froo + «)froo) | r=Q = , 

("I)*" I,-., 



(1.6) 



r - r 



It is given explicitly by the formulae, 

r + t- r 



(v £ _) hee (t,r) =P_ [r + t, 



r-t — ro 



P- t 



(ui) free (i,r) =P+ [r-t, 



The pulse (vi)f rC o corresponds to an incoming wave, and (v £ + )f rcc is the sum of 
two outgoing waves, one from P + , and the other from the focusing of the incoming 
wave. 
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Theorem 1.1 (Sub-critical case). Assume that a > max(0,p — 2). 

• If a > p — 2 > 0, then there exists eq > such that for any e G]0, eq], ( |l.5| ) /ias a 
unique, global, solution v £ G C ([0, oo[xR+). Moreover, the following asymptotics 
holds in L°°(M + xl+), 

v s ± (t,r) = (i4)froe(t,r) + O ( £ "««(i^+2-rt) , 
efl t t;|(t ) r) = ea t (4) frcc (t,r) + O (V^CW-p)) _ 

• If a > and 1 < p < 2, let T > 0. TTien i/iere exists e(T) > smc/i i/iai /or 
any e g]0, e(T)], (|l.5| ) /ias a unique solution v £ S C 1 ([0,T] x R+). Moreover, the 
following asymptotics holds in L°°([0,T] x 

4(t,r) = («|) free (t,r) + C»(e min ( 1 ^) (O (e + e Q | loge|) tfp = 2), 
e9 t u|(t,r)= £ a i ( 1 ;|) fr ee(t,r)+o( £ min ( 1 ^) (O (e + E a \ loge|) p = 2) . 

Theorem 1.2 (Super-critical case). Assume that < a < p — 2 or a = = p — 2. 

• If the equation is dissipative, a > 0, i/ie pulses are absorbed before reaching the 
focus. IfT>r a , 

limsup (||«i(T)|| £ -, ( o< P <T) + IKCT)||£-(0<r<T)) = 0. 

More precisely, for A > 7 define T(A, e) as follows; ifO < a < p — 2, then T(A, e) := 
ro — zq£ — \E a ^ p ~ 2 \ and if a = = p — 2, i/ien T(A,e) := ro — zo£ — A. For any 
T = T(e) > T(X,e), 

lim limsup (\\v e _{T)\\ L <*,( Q < r<T \ + ||v+(T)|| L =o( < r < T) ) = 0. 

• If the equation is accretive, a < 0, there exists T* < ro such that the family 
(v £ _,v £ + ) is not bounded in L°°([0,T*] x K + ) 2 . 

Remark. For < e < 1 and A positive, T(A,e) < ro — zo£, so the first part of 
the theorem shows that the absorption mechanism takes place before the incoming 
wave reaches the focus. Indeed, the pulses are initially supported in {|r — ro| < zqe}, 
so by finite speed of propagation, they do not reach the origin before t — r — z s 
(see Fig. |l|). Notice that in the dissipative case, our estimates are in {t > r}, which 

t A 



r - z £ 
T(X,e) 




Figure 1. Geometry of the propagation in the super-critical case. 



includes the focusing region, and its domain of influence. The key to the proof of 
Theorem 1.2 is to construct and approximate solution which is more accurate than 
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the approximation of nonlinear geometric optics and which permits us to penetrate 
with high accuracy to small distances from the focal point r = 0. 

Remark. If we considered wave trains, that is Uj(r, .) periodic instead of compactly 



supported, Theorem 1.1 would still hold. The proof we give works in both cases. On 
the other hand, the proof of Theorem |1.2| relies on the compact support assumption. 
We construct approximate solutions that solve ordinary differential equations along 
the rays of geometrical optics (see Sect. |^), and can be computed explicitly. In the 
dissipative case, these approximate solutions are absorbed before they stop being 
good approximations, proving thereby the absorption of the exact solution. In the 
case of wave trains, the computation of the counterpart of these approximations is 
a project for the future. 

Remark. As we mentioned above, in the supercritical cases, our framework is re- 
stricted. We assume tha t th e coupling constant a is real, while we do not make this 
assumption in Theorem 1.1. It would be interesting to know what happens when, 



for instance, a is pure imaginary; no absorption can happen, for the equation in that 
case is conservative. A partial answer is given in IKX04 in the case a = = p — 2, 
on a system which is a simplified model for (|0|): an arbitrary phase shift appears, 
varying like loge. In the supercritical framework, a — and p > 2, one may expect 
even more pathological behaviors. 

In Section |[ we prove two stability results. In Section ||, we discuss existence 
results, and establish estimates for the sub-critical case. Theorem 1.1 is proved in 
Section ||, and Theorem 1.2 is proved in Section g. 



The results of Theorem 1.2 were announced in [CROlf 



2. General stability results 

In this section, we state two general approximation arguments. The first one will 

and the second one, Theorem 1.2. Our first result 



allow us to prove Theorem 1.1 



is an easy estimate, proved in [CR04] 



Definition 2.1. For t > 0, we denote by r*i (resp. r+) the set of all speed minus 
one (resp. plus one) characteristics connecting points on the initial line {t = 0} to 
points at time t. We also denote T* = T*_ U r + . 

For a characteristic 7 £ T*, we use the convention that 



/ 



stands for the integral of / along 7, parameterized by the time variable, 



/= / f(s,r(s))ds, 



where (s, r(s)) sE [o,t] is a parametrization of 7. In particular, if / is nonnegative, 
then so is the above integral. 



Lemma 2.2 ([CR04|, Lemma 2.1). Suppose that w and f = (./+,/-) are bounded 
continuous functions on [0,T] x [0, oof satisfying in the sense of distributions 

(d t ± d r )w± = f± , w+{t,0)+w-{t,0) = 0, forO<t<T. 

Denote by 

M ± (t) := || W± (i)|U~ ([0 ,oo[). 
Then for < t < T one has 



M±(t) < max{M+(0),M_(0)} 



sup 

7-er'_ 



I/- 1 



sup 



1/4 
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To prove Theorem 1.2, we use a result in the spirit of GronwalPs lemma. The 
assumption < T < 5 guarantees that the support of the solution does not touch 
{r = 0}. 



Proposition 2.3. Suppose thatO<T<S, and w = (w- , w + ) G C H L°°([0, T] x 
R+) satisfies 

(d t ± d r )w± = f±(t, r)(w- +W+) + S ± , 
w±\t=o = u> ±, 
with supp wq± c [5, +oo[. Suppose that 



(2.1) 



(2.2) 

Then 



C l = C l {f) 



/ sup \f_\dt+ / sup \f+\dt < oo. 
Jo 7-er*_ Jo 7+er* 



sup ||w± 



0<t<T 



\u>o±\\l<~ + sup 



7±ei± J~f± 



\S±\ 



with C 2 = max(Cie 2Cl ,C 1 2 e 3Cl )- 



Warning. Note that in hypothesis ( |2.2| ), the supremum is inside the integral. The 
estimate would not be true with the supremum outside. 

Proof. Let (i, r) G [0, T] x]0, oo[, and denote 7- = 7-(i, r) the characteristic from 
(0, i + 7') to (i, 7'). Duhamel's principle for W— reads 

w-(t,r) =w -(t + r) + f-x(w-+w+)+ S— 

J "f _ t/ 'V— 



Gronwall lemma, along with assumption (2.2), yields, for any t G [0, T], 
||«>_(t)||Lo= <e Cl [ + sup 



(2.3) 



sup / \f-.w + \ + sup / \S-\) 



<e Gl ||tuo_|| L oo +Ci sup ||u; + (a)||L» + sup / |S_ 

V 0<s<t 7-Gr*_ ^7_ 



Similarly, 

(2.4) |K(f)|| £ « <e Cl + sup / sup / |.S' + | 



7+ er^ J 7+ 



7+ er^ J 7+ 



For the w_ integrals on the right in (2.4), use estimate (2.3) to find 



\f±w-\ < |/±| e Gl ||«;_(0)||+Ci max \\w+(r) 

\ 0<T<S 



where j± is parameterized by the time s G [0,t]. Introduce 

m±(f) := sup ||w±(s)|| , 

0<s<t 



sup 

7- erf. J 7 _ 



and take the supremum on r in (2.4) to obtain 



K(*)ll < ^ie Cl (E(ll w ±(°)ll+ SU P / l 5 ±l) + 

V ± V 7±er' ± J 7± J 

+ sup / |/+(s,r)| m + (s)ds + sup / |/-(s,r)| m+(s)ds . 
7+ er*.^7+ 7-eriJ 7 _ / 
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Therefore, 



|m + (t)|| < C ie Cl (^(|| W± (0)||+ sup / |5±|) + 

V ± V 7±er* ± J 7± / 

+ sup / \f + (s,r)\m+(s)ds+ sup / \f-(s,r)\ m+(s)ds 
7+er* + J 7+ 7 _eriJ 7 _ J 

^ Cie Cl ( W|k±(0)|| + sup / |5±|) + 

V ± V 7±er t ± J 7± / 

+ sup \f + (s,r)\ m + (s)ds + sup |/_(s, r)\ m + (s)ds . 
J a 7+er* + Jo 7_er*_ y 



Applying Gronwall's lemma, using assumption (2.2), yields 



m+(T) < de 2Cl 



K(0)|| + ||™_(0)||+ sup 

7-er'i' J~i 



\S- 



sup 

7+er^ J 7+ 



|S 4 



This inequality, along with ( |2.3| ), proves the proposition 



□ 



The following example contains the core of the proof of Th. 1.2 
Exam-pie 2.4. Consider the case 

f e ± (t,r) =s a r 1 -P\vl(t,r)\ p - 1 , 
where v e _ solves an initial value problem of the form (compare with ( |5.lD below), 

(d t - d r )W_ = F e (51) , trl |t=0 = P_ (r, ^ 

Recall that suppP_(r, .) C [— z , zq]. For Wq ± , consider pulses with the same sup- 
port as ufi t _Q. Then by finite speed of propagation, we can take 8 e = r — z e in 
Prop. 2.2 . For < t < 8 s , the maximum of /|. at time t is estimated by 

F a r 1 ~ p \\v 6 ll p_1 

fc L ll w -lli«>(o<t<t)j 

where r is such that r + 1 = 6 e = ro — zq£. Therefore we have 

sup \fl\dt+ f~ sup \fl\dt<2\\¥\\ p -J (Q<t<t) -i) 1_p <tt 
cr* .Iq - y+ er t Jn 



o 7_er'_ 



< C p e a \\v'L\\ p L J( < t <t- j x 



(S e -t) 2 - p ifp>2, 
\log(S s -t)\ i{p = 2. 



Thus, when v e _ remains bounded and T e is chosen so that e a (5 £ — T e ) 2 ~' p (resp. 
e a \ log(<5 £ — T e )\ if p = 2) is bounded independent of e, then we can use Prop. pT^ . 
This is the case in particular if T £ = T\ e defined in Th. 1.2. 

3. Existence results 

In this section, we prove two kinds of results concerning the existence of solutions 
to (1.5); local existence in L°° before the wave meets the boundary {r = 0}, and 
global existence in W 1 ' 00 when the boundary condition has to be taken into account. 
The reason appears in Th. [hi] and 1.2; in Th. 1.2, the phenomena we want to prove 
occur before the pulses reach the boundary, while Th. 1.1 includes the caustic 
crossing. 
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In the first case, we are interested in a problem 

j (d t ±d r )v e ± = e a r 1 - p g(v e _+v s + ), r>0, 

where suppug T C [<5 e ,+oo[ for some 5 s > 0, and Vq t G L°°(R+). Then by finite 
speed of propagation, the term r 1— p is harmless at least up to time 8 £ . In that case, 
local in time existence of solutions to is easy. 

Lemma 3.1. Fix a>0,p> 1. Let Vq t G L°°(K+) such that suppvp-p C [<5 e ,+oo[ 
for some d E > 0. Then there exists T £ , with < T e < 8 e , and a unique solution 
{v e _,V%) G L°° n C([0,T e ] x M+) 2 to the initial value problem 

When the incoming wave v £ _ reaches the origin {r = 0}, a boundary conditi on is 
needed in order to solve the above system. We are interested in that given in ([O]), 
that is 

(«- + <)U = °- 

Consider the mixed problem 

r (d t ± d r )v e ± = e a r l - p g{v e _ + v e + ), r > 0, 
(3-2) | (fi+<)U = 0, 

The boundary condition compensates the singularity r x ~ p when r gets close to zero. 
Indeed, Taylor's formula yields, for C 1 solutions, 

(vt + v e + )(t, r) = rd r (v e _ + v e + )(t, r) + o(r), as r -> 0. 

Now from the differential equation, we also have 

d r (v s _+v e + ) = d t (v s _-v s + ), 

so if we know that the time derivatives of v £ remain bounded, then the singularity 
r 1_p is compensated. We have precisely, 

Ar 

(3.3) \(vt+v e + )(t,r)\ <__(|t,i| + | v «| + | e ^t;L| + | e ^|)(t,r), Vi,r>0. 



This is the strategy we used in [CR04], Proposition 3.4, to prove local existence, in 
the case p > 2. Notice that at this stage, the dependence upon e is unimportant, 
and that p > 1 suffices for local existence. 

Lemma 3.2. Fix a > 0, p > 1. Let v^ T G Vr ll0o (M+) such that 

(3.4) v e „ (0) + v e Q+ (0) = , and d r v e Q _ (0) - d r v e 0+ (0) = 0, 

then there exists T e > and a unique solution [v e _,v%) G C 1 fW^QO, T 6 } x R + ) 2 
of ©■ 

As recalled in the beginning of this section, such a result will be needed only 



in the proof of Theorem 1.1, and not in the proof of Theorem 1.2. From now on 



in this section, we assume a > max(0,p — 2). In |CR04|, we also proved that the 
solutions of (3.2) with e = 1 and p > 2 are global provided that the initial data 
vq± are sufficiently small. 

Proposition 3.3 ( ]CR04[ |, Proposition 3.5). Fix a > 0, p > 2. There are constants 
K\ and K[ > so that for all initial data ipo G C 1 ([0, oo)) satisfying 

||^o,<9,.-0o||l~([o,oo[) < K x , 

and the compatibility conditions 

^o+(0) + ^o-(0) = 0, and 9 r ^ +(0) - 3 r Vo-(0) = 0, 
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there is a unique solution ip G C 1 ([— 00,00] x [0, 00 [) of 

((d t ±d r )^ ± = r 1 -f 5 (^_+^ + ), 

(3.5) J V-M) + V+(*,o) = 0, 

In addition, 

°°([- 00,00] X 

The idea is then to find a scaling such that we can use the above proposition to 
prove global existence for (1.5) when a > p — 2 > 0, along with useful estimates. 
Try 

(3.6) v e ± (t,r)=e^l(r,p)\tr , 4> e ± (t, p) = £^4(er, ep) . 

Then ip± solve the differential equations, 

(d T ± d P W± = e Q (ep) 1 - p e (p - lh <?(C + 4>+)- 
Choose 7 so that the powers of e cancel, 

a+1— 1)7 = 0^7 = 1 — - . 

P - 1 

Since 7 is negative, 

||^|(o)|U~ , ||^|(o)|| £ -— *o, 

so we can apply Prop. 13.31 , for e sufficiently small. Moreover, it provides L°° es- 
timates for We deduce that d r v± e L°° from the differential equations 
and (3.3). 

Corollary 3.4. Assume a > p — 2 > 0. Then there exists Eq > suc/i £/iai /or 
< e < e , (1.5) /las a unique solution (vt,v+) 6 C 1 nW 1 '°°(R + xM + ) 2 . Moreover, 
there exists C such that for any e £]0,£o]> 

||«±,£d t w!l| i0 o (R+xR+ ) < C. 

We now have to prove global existence when 1 < p < 2 and a > 0. Since 
local existence is known (Lemma 3.2), global existence is a consequence of a priori 
estimates, which follow from (3.3) and Lemma |2~2] . Define 

m|(i) = sup (||*4(s)|U» + \\ed t vUs)\\ L ~) . 

0<s<i 



By assumption, m§_(0) are bounded independent of e €]0, 1]. Lemma 2.2 and ( J3.3| ) 
yield, for T > 0, 

ml(T) + m^(T) < C (m £ _(0) + m* (0) + h p (e,T) (m £ _(T) + m^(T)) p ) 
< Co + C (h p (e, T) (mi(T) + m £ + (T)) p ) , 

with 

£ a T 2-p ; if 1 < p < 2 , 

e «l og (l + Z) ,ifp = 2. 

In particular, for fixed T > 0, h p (e,T) — > as £ — » 0. Therefore, for any T > 0, 
there exists e(T) > such that for any e e]0, e(T)], w e exists and, 

m e _(T) + xn E + (T) < 4C . 
If mI(T) + m^(T) < 4C , it follows that 

mi(T) + m £ + (T) < C + CTi p (e, T)(4C ) P , 
and therefore for < e < e(T), 

ml(T) + m £ + (T) < 2C 



10 



R. CARLES AND J. RAUCH 



This proves that in fact mf_(T) + m+(T) < 4Co since if that were not true, there 
would be a first T* where mi (T) + m £ + (T) = 4Co, and at that value of T the above 
estimate leads to the contradiction ACq < 2Cq. 

Proposition 3.5. Assume that a > and 1 < p < 2. Let T > 0. TTien there 
exists e(T) > such that for < e < s(T), ( |1.5| ) /ias a unique solution 

(vt,v £ + ) e c 1 nf^do.T] x M + ) 2 . 

Moreover, there exists C such that for any e €]0, £(T)], 

||v|,eQ{t!|.||i,«.([o ) T]xR + ) < C< 



4. The subcritical case 



Now Theorem 1.1 is a straightforward consequence of Lemma 2.2, Corollary 3.4 
and Proposition 3.5. In the statement of Theorem we distinguished three cases; 
p > 2, p = 2 and 1 < p < 2. The distinction p<2 appeared in the previous sec- 
tion, in Corollary \A and Proposition |3.5| . It corresponds to the question of the 
integrability at infinity of the mapping r i— > r x ~ v . The further distinction p = 2 
corresponds to the local integrability of this mapping. 



(4.1) 



Define the remainder w± = vfj_ — (u±)free- It solves the mixed problem, 
f (d t ±«9 r )u4 = £ a r l - p g{v e _ + v s + ), 

(™L + <)| r=0 = o, 

el -\- T> ( r ~ r ° " 



Let / be an interval of the form [0, T[, with T £ R + U {+oo}. From Lemma 2.2 
have 

||«4IU-(/xr + ) <2e||Pi|| i oc+2supsup / eV~P ff («I+^). 



we 



Using (3.3), we also have, 

(4.2) |K ± |U~ (/xR+) < Ce + C (sup sup / /° ) K'^IIWx**) 



Differentiating (4.1) with respect to time yields, using the differential equation (4.1) 
to find the initial data, 

( (d t ± d r )e.dw £ ± = e a r l -Pg'{vt + v s + )ed t {v E _ + v%), 

:d t (w s _+wl)\ r=0 = 0, 



edtw 



T lt=0 



eiedrP! + d z P!) r, 



Since the initial data for v± are supported in \r — ro | < zo£, the term r 1_p in the 
initial data for edtw^ is harmless regarding to L°° estimates. From Lemma 2.2 and 
(I3.3D, we have, 



|]^<4IU~(/x R+ ) < e^ 1 -") + sup sup / ll«Vfif>l£« (IxR+) 

\ te/ 7 er' J 7 l r + e ) p I 



7 ei 1 j 7 

The case p > 2. Assume a > p — 2 > 0. Then we have 



sup sup . 

t>o 7 er* J 7 (»" + e) p 1 



< Ce 



Q+2-p 
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and from Corollary 3^, there exists eo such that if e G]0, £o]j ||^ E ; £d t v E \\ 1 
is bounded. This proves the first part of Theorem 



1.1 



The case p = 2. When a > p — 2 — 0, then from Prop. 3.5, for every T > 



there exists s(T) such that if e e]0, e(T)], ||w e , edtV e \\ La 
Moreover, for any fixed T > 0, we have 



'([0,T[xR+) 



remains bounded. 



sup 



7 



— < Ce a log - . 

■ e e 



The case 1 < p < 2. When a > and 1 < p < 2, then the only difference with the 



previous case is that the mapping r 

e° 

sup 

7 er T 



< sup 



is locally integrable, hence the bound 



, 7 (r + e)P- 1 
This completes the proof of Theorem |l.l[ 



< Ce a T 2 - p . 



□ 



5. The supercritical case 



We conclude by proving Theorem 1.2, using Proposition 2.3 



Before going into details, we explain how to construct the approximate solutions 
that lead to the result. Assume for instance that a > and, since we are in a 
supercritical case, p — 2 > a. The hypothesis a > suggests that the nonlinear 
term in ( |1.5 ) is negligible when r is not too small. This is the outline of the proof 
of Theorem 1.1, and we could prove this way that («j_)f ree is a good approximation 
for v± at least for tq — t 3> e~f^ . This boundary layer is larger than in the critical 
case a = p — 2 > studied in [CR04], and nonlinear effects possibly occur sooner. 



Considering the case a = gives us a further hint. We proved in |CR02] that a 
good approximation at leading order was given by 



npp) 



that is, solutions of ordinary differential equations along the rays of geometrical 
optics. A natural generalization of this approach to the case a > leads to 



(5.1) 



(dt±d r )(v%) app =eV 

(U±)app| t=0 = P± 



- p g((v e ± ) app ), (t, r) G ([0, r - z e[xR* + ) 
ro ' 



We consider the region t < ro — zqE so that no boundary condition is needed on 
{r = 0}, since the compact support of P± and the finite speed of propagation make 
(vfjjapp zero in the region {r + t < r — z e} (see Figure l]). 
Introduce the function F p , defined for y > x > by, 

L -<x 2 -p -y 2 -?), ifp>2, 



F P {x,y) 



ds 

s'P- 1 



P 



II 



Then (5.1) can be solved explicitly, 



(ul)app(^r-) = 



(«+) apP (*,r) = 



log-, lip 

x 



P-(r + t,z) 



(l + a2-P(p - l)s a F p {r, r + t)|P_(r + t, z)^- 1 
P+(r-t,z) 



(l + a2-p{p - l)e a F p (r - t,r)\P + (r - t, z)^- 1 ) : 



_ r+t-rp 



_ r-t-rp 
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We have, 

supp(wl) app = {\r + t - r \ < z Q e}, supp(i^) app = {\r-t-r \ < z e}. 
Since for any y > and any p > 2, F p (x, y) — ► +oo, the sign of a is crucial. 

• If a > (dissipative case), then (i>i) a pp tends to zero before reaching the 
focus. 

• If a < (accretive case), then both (i>i) app and (i/j_) app blow up in finite 
time. 

In the accretive case, (w^_) app may blow up at time T* < ro, that is, before any 
focusing. To avoid useless distinctions, we assume P+ = 0, so that the only "inter- 
esting" phenomena occur when t approaches ro. 

The explicit formulae show that when p > 2, for s a r 2 ^ p 3> 1, 

(4) app (i,r) = (4) frcc (i,r) +0(e a r 2 -P). 

Assume first p — 2 > a > 0. Let A > 0. We prove that (t>§.) app gives a good 
approximation of the exact solution in the region r > Ae^ 5 before the focus. 
Notice that in this region, (wj_) app and (v±)f rcc have ceased to be close to each 
other; the nonlinear effects are significant. By finite propagation speed, this area 
is defined by 

t < r - Xe~ - z Q s = T\ jS . 
We prove that (wj_) app remains a good approximation of the exact solution up to 
time Tx, e . 

Proposition 5.1. Assume p — 2 > a > 0, and let A > 0. Then (v±) app is a good 
approximation of the exact solution at least for t £ [0,T,\ jE ], 

||«±-(«±).p P |U~([o,T*..]xB + ) = 0(e 1- ^) (0(s) ifp = 2). 



Proof of Proposition 5.1. The proof relies on Proposition |2.3[ which we apply as in 



Example 2.4 



By finite propagation speed, both the exact and approximate solution are zero 
in the region 

{(t,r) £ %T x ,e] xK+; r + t < T A , e }. 
Define the remainder w± — v± — (wj_) app . Before rays reach the boundary {r = 0}, 
( 9t ± d r )wl = eV"* (g{vt + v%) - 9((«l)ap P )) , 

w ±| t=0 = T£Pi (r,^ 
In order to apply Proposition |2.3| , write the right hand side as 
eV"* (g(vl + <) - .g((4)ap P )) =e a r 1 ~P {g{vl_ + v%) - ) app + «) app ) 

+ S(0-)app + «)app) ~ .g(M) ap p)) . 

Using Taylor's Theorem, the first term satisfies 

eV"* (g(vi + <) - g((vl) app + «) app )) = e V"*(ti£ + <)/ £ (i, r), 

where f E is uniformly bounded on any set on which the families v E and (u e ) app 
are uniformly bounded. The point now is that this term has exactly the properties 
mentioned in Example 2.4; in particular, the assumptions of Proposition [2.3| are 
satisfied up to time T = T\^ E . 
Define the source term 

S|(t,r) := eV-"( ff ((«I) am , + «) ap p)- ff ((«i) am ,)). 
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Then supp<S± C supp(w^-) ap p , and 

' {d t ±d r )w% = e a r 1 - p f e {t,r){w e _ + w e + ) + S £ ± , 

r - r s 



(5.2) 



w ±\t=o = ^ ePl I r ' 



From the explicit expression of (vfj_) app , we see that for any fixed A > 0, there exists 
C\ > independent of e such that 



(5.3) 



04) 



±; a PP|lL°°([0,7\, e ]xR + ) 



< c x . 



We shall prove that there is e(A) > 0, so that for e £]0, e(A)] 



(5.4) 



|W ± ||Loo([ 0jTA _ e ] xR+ ) 



< Cxe 1 - 



(C x e ifp = 2 



This implies the error estimate of Proposition 5.1. For the sake of readability, we 
will omit the distinction p = 2, and keep the notation e ~ f~ 2 , with an obvious 
convention. 

From Lemma 3.1, v £ , hence w £ , is defined, bounded and continuous, locally in 
time. At time t = 0, it is of order e, so there exists t £ > such that 

||w±l|i«([o ) t«]xK + ) < 2e||Pi||ioo, 
and we have, possibly increasing the value of C\, 



P±IU~([0,t=]xK + ) 



So long as w £ ± is pointwise bounded by 2C\, f £ remains uniformly bounded. As 
noticed in Example 2.4, there exists (72(A) such that 

"" supeV-* |/ £ |dt < (7 2 (A). 

Now consider the source terms. On supp(u!j_) app (of size 2zo£, and transverse to 
any 7 6 e '), r > ro/2, therefore the singular term r 1_p is bounded and 



sup 

7 _ er T(A, E ) J 1 



\S £ _\ <Ce 



l + Q 



More delicate is the treatment of S £ + . Because supp(ul) app is of size 2z e and 
transverse to any 7+ £ rT , one has 



sup 
If p> 2, 



\S £ + \ < Ce a 



\si>- 2 +2z e 



Cs a F p (Ae^ , Ae^ + 2z ej 



(p-2)e a F p (\e^,Xs^ + 2z e) =e a MAsAj P - (Ae^ +2z e 



2-p 



2-p\ 



The case p = 2 yields the same estimate. From Proposition 2.3, u>|_ remains 

pointwise bounded by 2(7 A for f £ [0,T(A,e)], for e £]0,e(A)]. This yields 

and completes the proof of Proposition |5.l[ □ 
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Proposition 5.1 describes the behavior of the exact solution v £ up to time T(X, e). 
The outline of the end of the proof in the dissipative case is as follows. At time 
t = T(\,e), the approximate solution (v £ ) app is of order \(p~ 2 ) / (p^ 1 ) if p > 2 
(1/| log A | if p = 2). Letting e go to zero, with A sufficiently small, shows that v £ 
becomes arbitrarily small when approaching the focal point. Since the equation is 
dissipative, this means that v £ is absorbed. The end of the proof of Theorem |l.2| 
relies on energy estimates. For q > 1, define 



q\s\ q 1 sgns. 



Then g q is a non-increasing odd function of s which is homogeneous of degree q. 

In ( |l.5| ), multiply the equation satisfied by v £ _ by <? g _i(ui), and the equation 
satisfied by v £ + by g q -\(v £ + ). Summing yields 

d t {\vt\ q + \v £ + \ q ) + d r (\v £ + \ q - \v e _\ q ) = 

- a2~Pr 1 -P + <fe-i«)) 9 P (vL + v £ + ) / {p + 1) . 

The signs of both g q -i(v £ _) +g q -i(v s _) and g p {v £ _ +v £ + ) are equal to the sign of the 
larger of Therefore, when a > (dissipative case), 



(5.5) 



dt(\vt\ q + \v £ + \ q ) - d r (\v £ _\ q - \v £ + \ q ) < 0. 



For a fixed t > 0, integrate this inequality from r = to r — t. Recall that when 
r = 0, we have \v £ _ \ = \v+\, so this yields 

d t (\v e _\ q + \v £ + \ q )(t,r)dr - \v £ _\ q (t,t) + \v £ + \ q (t,t) < 0. 



o 



Therefore, 
(5.6) 



dt I (\v £ _\ q + \v £ + \ q )(t,r)dr <2\v £ „\ q (t,t) 



Let A > 0. From Proposition 5.1, for < t < T x , e , 

\vl\ q (t,t) = \(v e _) app \ q (t,t)+0{e< 1 - 

For t > ^±^, (vt) app (t,t) = and 

\vt\ q (t,t) = O (e*( 1_ 3£s)) . 

Let T > T A , £ and integrate ( |5l| between t = T A , e and i = T. 

ll^(^)ll^(o,T) + \K(T)\\Uo, T) <\\v S -(TxMl H0 ,T x , e ) + \K( T ^)\\Ua, T ^ 

+ C(\) q Te q ( 1 -^\ 

and 

\\vt(T)\\ L <,(p t T) + IK(T)|| L , (0 ,T) <\\vl(T X ,e)\\ L i(0,T x , e ) + ||«+(TA, e )|U*(0,T\, s ) 

+ C(A)T 1 / 9 e 1 -^. 

Letting q — > oo yields 

IKCr)IU~ (0 ,T) + KCOHioofcT) <ll«I(Tv)IU-(o,7\,.) + IK(7A, e )IU-. {0 ,^,.) 

+ C(A)e 1 "^. 



Using Proposition 5.1 again, and the fact that supp(w^_) ap p C {\r — t — ro\ < zqe}, 
we also have 

ll«ICniU«<0,T) + \K(T)h~(.0,T) < ||(«i)app(T A , e )|U»(0,T,, s ) + C(A) £ 1 -^. 
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From the explicit expression for (vL) app , 
(5.7) W(vlU P (Tx,e)h~(o,T x>e ) <Cx ' 



Xp- 1 , if p > 2, 
1/| log A| , ifp = 2, 



where C does not depend on A. Thus if p > 2, 

IKCOIU-flvr) + IK(T)|| £ » (0i t) < CA^? + C(A)e 1_ '^2 . 
Therefore, for any T > ro, and any A > 0, 



(5.8) limsup (\\vl(T)\\ Lao(0tT) + |K(T)||l~(o,t) < CX~ 



Letting A — > yields the first part of Theorem 1.2 for p > 2. The case p — 2 is 
straightforward. 

When the equation is accretive (a < 0), the energy estimate ( |5.5D becomes 

(5.9) fcfluil* + \v%n(t,r) - d r (\vt\ q ~ \v s + m,r) > 0. 

Since we assumed P + = 0, (ul) app blows up in finite time, while (i>5_) app does not. 
The mechanism occurs quite the same way as the cancellation of (wi) app in the 
dissipative case. For a fixed t > 0, integrating (^^) between r = and r = oo 
yields, by finite speed of propagation, 

/>oo 

d t (\vl\i + \v^)(t,r)dr = d t / (|»L|9 + K.|9)(i,r)dr>0. 

Jo 

Let ro > T > Ta,s and integrate the above inequality between £ = Ta, £ and t = T . 

Wv-{T)\\l m+) + \K(T)\\l m+) > \\vl(T x , E )\\l q(R+) + K(T Aie )||« s(B+) . 
Letting q — > cxo yields, 

IKCTJIUoo^) + ||<(T)|| L » (R+) > ||t,i(T A:e )|| L ^ (R+) . 
From Proposition 5.1, it follows 

Uminf (lluiCT)!!^^) + \\v e + (T)\\ L ^ {m+) ) > liminf ]](«i) a pp(TA, e )]]i,oo (H+) . 

Letting A — > yields the last part of Theorem 1.2 with T* — ro, using the explicit 
form of (i>i) app . As we already mentioned, the result may hold with T* < ro if 
(f4-) ap p blows up before (fi) app , in which case the proof is essentially as above. 
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